Locally extracting scalar, vector and tensor modes in cosmological perturbation theory 
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Cosmological perturbation theory relies on the decomposition of perturbations into so-called 
scalar, vector and tensor modes. This decomposition is non-local and depends on unknowable 
boundary conditions. The non-locality is particularly important at second- and higher-order be- 
cause perturbative modes are sourced by products of lower-oder modes, which must be integrated 
over all space in order to isolate each mode. However, given a trace-free rank-2 tensor, a locally 
defined scalar mode may be trivially derived by taking two divergences, which knocks out the vector 
and tensor degrees of freedom. A similar local differential operation will return a pure vector mode. 
This means that scalar and vector degrees of freedom have local descriptions. The corresponding 
local extraction of the tensor mode is unknown however. We give it here. The operators we de- 
fine are useful for defining gauge-invariant quantities at second-order. We perform much of our 
analysis using an index-free 'vector-calculus' approach which makes manipulating tensor equations 
considerably simpler. 



I. INTRODUCTION 



Perturbation theory in cosmology rests upon the decomposition of perturbations into scalar, vector and tensor 
parts, on a background which has constant curvature. This is a generalization of Helmholtz's theorem to tensors on 
3-spaces of constant curvature [1-3] . This split is usually performed non- locally: either in harmonic space or using an 
implicit integral over a Green's function (they are equivalent). Because of this, we must specify boundary conditions 
■ in order to define these non-local variables. Furthermore, we must know a perturbation variable everywhere in order 
to specify a given scalar vector or tensor type of perturbation somewhere. For example, let us consider Hclmhotz's 
theorem in 3-d fiat space. Any vector V can be trivially written in terms of a scalar ip and vector A: 



> 

in 

Tj- ! V = Vip+cur\A, (1) 



where 

V 2 V = divF, (2) 
V 2 A = -curlV, (3) 

which follow from the standard vector calculus identities div curl V = and curl curly = — V 2 y + Vdivy (in 
Euclidean space). Unique solutions for ip and A exist provided they vanish sufficiently fast at infinity. That is, ip 
and A are inherently non-local, requiring knowledge of V everywhere just to be defined at a single point. On the 
other hand, we can think of <p — div V as a pure scalar degree of freedom which is defined locally wherever V is; 
similarly, curl V is a pure vector degree of freedom. So, given V we can isolate locally defined, unique, scalar and 
vector degrees of freedom by differentiating it appropriately. This is useful in electromagnetism where, for example, 
it is common to write down decoupled wave equations for the scalar and vector potentials which arc equivalent to 
Maxwell's equations. Of course, the solution to these equations involves the same reliance on boundary conditions 
as defining ip and A do; yet the split of the dynamical wave equations themselves into decoupled solcnoidal and 
irrotational parts is conceptually useful - and this doesn't have to rely on non-local conditions, or integrals over all 
space. 

Given the unique nature of cosmology, it seems perverse for us to only be able to define gravitational waves or the 
gravitational potential utilising unknowable boundary conditions and conditions outside our horizon. The universe 
could be spatially infinite, implying that we cannot arbitrarily assign boundary conditions at infinity (which could 
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depend on the start of inflation, for example). For scalar modes, however, we can reformulate perturbation theory 
into an equivalent local theory because we can simply take divergences to form objects which are locally defined pure 
scalars. Because of this they have a meaning which is well defined physically. Similarly for vectors. But what is 
the equivalent operation for tensor modes? Given a rank-2 tensor, how do we locally isolate the tensor degrees of 
freedom? That is the goal of this paper. For readers who just can't wait to find out, the answer is given by Eq. (37). 



Notation and identities 



We assume an FLRW geometry of curvature K, and all relations are defined on this background. Under perturbation 
at order n, all relations below hold for objects of perturbative order n; for lower perturbative order, the commutation 
relations below have curvature correction terms added to them. Furthermore, we assume the objects we deal with are 
tensors in the full spacetime, which implies that they must be gauge-invariant objects at order n from the Stewart- 
Walker lemma [f]. 

Given the usual 4- velocity u a wc define the spatial metric h ab = g ab + u a u b and volume element e abc = u d rj abc d [5]. 
All other rank-f and -2 tensors used here arc orthogonal to u a , and rank-2 tensors are projected symmetric and 
trace-free (PSTF) which we denote using angle brackets on indices. 1 We define the conformal (comoving) spatial 
covariant derivative acting on scalars or spatial tensors as V a = ah b V b — aD a , where a is the scale factor and D a is 
the spatial derivative normally used in the covariant approach. 2 We use V a as it commutes with w a V a and is the 
covariant derivative normally used in the metric approach to perturbation theory when its index is downstairs (when 
K = it is just the partial derivative); all expressions below are covariant, and indices are raised and lowered with 
g ab (or h ab for PSTF objects). The irreducible parts of the spatial derivative of PSTF tensors are the divergence, 
curl, and distortion, defined as [4] 

(4) 
(5) 
(6) 

Then, the spatial derivative of a rank-n PSTF tensor X An = X aia2 ._ an may be decomposed as (for n = 1, 2, 3) 3 



divA\.. c 




curl X ab ... c 




d\sX ca ,,, b 





2n — I n 
V b X An = d\M X {An l h an)b — cur\X c{An l e an)b c + d\sX bAn . (7) 



Note that the divergence decreases the rank of the tensor by one, the curl preserves it, while the distortion increases 
it by one (and all are PSTF). Keeping this in mind one can drop the indices on differential operators as long as it's 
explicit the valance of the PSTF tensor which is being acted on (though this can sometimes be confusing). This 
considerably simplifies the appearance of the equations. 

There are important commutation relations for each of the three invariant derivatives which we need, which are 
derived from the 3-Ricci identity. The Riemann tensor on a 3-space of constant curvature is given by: 

i OK 

3 R ab cd =- 3 Rh a [c h b d] = -- ¥ h a [c h b d] => V [a V b] X c = KX [a h b]c (8) 



1 Because of the confusing lexicon of scalar, vector and tensor modes, we will try to call things like a 'vector' (a spatially projected 
4-vector, which in an FLRW background is a 3-vcctor) a rank-1 tensor, and a 'tensor' a 'rank-2 tensor', etc. Ridiculous, but it should 
avoid confusion. 

2 Strictly speaking D a Xi, c = h a a h b b ■ ■ ■ h c c V a /X(,> ... c /, which we must use for objects which are not at maximal perturbative order. 

3 All our derivatives may be generalised to projected covariant derivatives in a general spacetime by replacing V 4 D. Then we have 
the invariant parts of derivatives of PSTF tensors div , curl & dis which we discuss in the appendix. These don't commute with the 
derivative along u a , even in an FLRW geometry. 
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Then, when acting on a rank-1 spatial tensor (a 3- vector) the following relations hold: 



div curl 




( y J 


dis curl 


= 2curl dis 


(10) 


rnrl 2 -i- V 2 


— VHiv 4- OK 

V Ul V \^ £iL\. 


(~\ 1 "> 


div dis 


= + S^div + K 


(12) 


V 2 div 


= div V 2 - 2Adiv 


(13) 


V 2 curl 


= curlV 2 


(14) 


V 2 dis 


= disV 2 +4/\dis 


(15) 



Here we have used the notation whereby V with no index represents the gradient of a scalar. These are the familiar vec- 
tor calculus identities, on spaces of constant curvature. To be clear about notation, the first of these is e abc V a W b X c = 
0, and the second V {a {s b)cd V c X d ) = e de{a V d V b) X e + e de{a V d V e X b) - |e d(o6) V d V f X f = e de{a V d W b) X e . Note also 
that the first is a rank-0 tensor equation (curl preserves rank, div reduces by one), while the second is rank- 2 (dis 
increases by one). 

Acting on any PSTF rank-2 tensor, the following commutation relations hold: 



curl div = 2div curl (16) 

curl 2 + V 2 = |dis div + 3iT (17) 

dis div = | div dis - |V 2 -5K (18) 

(disdiv)curl = curl (dis div) (19) 

div curl div = (20) 

V 2 curl = curlV 2 (21) 

V 2 div = divV 2 -4Adiv (22) 

V 2 divdiv = div div V 2 - 6/vdiv div (23) 



where Eq. (18) is given in [4]. Some of these are really tedious to derive, but are easy to use in this form. Note that 
both V 2 and curl 2 can be written in terms of dis div and div dis . 

II. THE DECOMPOSITION THEOREMS 

We discuss now the decomposition of rank-1 and -2 (PSTF) tensors into their SVT parts. In particular, we shall 
separate local SVT variables from non-local ones. We shall denote local SVT variables by a breve Our notation will 
be such that the same letter with a different number of indices represent different but usually related objects (e.g., 
S a = VaS* and S a b = V< 5 6 ) = V( a V b )5 etc.). 

A. rank-1 tensors 

For rank-1 tensors the scalar and vector parts correspond to the curl- and divergence- free parts respectively. Given 
a spatial rank-1 tensor on an FLRW background at maximum perturbative order, X a — S a + V a = V a S + V a , it is 
easy to form rank-1 tensor quantities which are pure scalars or vectors by the following rules: 

Scalar: S a = V a div X = V a div S = V a V 2 5 
= (V 2 - 2A') V a 5 

Vector: V a = curl X a = curl V a . (24) 

By taking spatial derivatives we have quantities which are scalars and vectors but which remain local. We can then 
formulate the non-local variables from the original tensor through these variables by the formal solution 

S = V- 2 divV = V- 2 V- 2 div£, (25) 

V a = - (v 2 - 2A') 1 curl 2 X a = - (v 2 - 2Kj 1 curl V a . (26) 
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Here we have used standard inverse Laplacian notation whereby ^V 2 — nKJ stands for the solution of the corre- 
sponding elliptic differential equation: assumptions about behaviour at infinity or boundary conditions must be made. 
While X a , S a and V a have compact support, S and hence V a do not [2]. 

Part of the conceptual utility of defining local scalar and vector quantities via Eqs. (24), which are the same rank as 
X a , is that the differential operations involved commute with the Laplacian V 2 , and time derivative. Therefore, if X a 
satisfies a wave equation with source, £[X a ] = S a , where C is a linear differential operator containing the Laplacian, 
then we may locally extract the covariant scalar and vector parts to find C[S a ] = V a div5 and C[V a ] = curl S a . 

One may further relate the local and non-local decompositions in Fourier space (which is inherently non-local, and 
relies on the functions having compact support - i.e., vanishing sufficiently fast at infinity). Defining a scalar harmonic 

basis in the usual way, V 2 q(S) = -k 2 Q<* s \ we find = -k 2 S<- k \ Similarly for vectors, V 2 Q a V) = -k 2 Q a V \ where 
we have two parities of orthogonal harmonics, (k 2 + 2K) 1 / 2 Q^ = curl <=> (k 2 + 2K) 1 / 2 Q ( - V) = curl Q^ v \ in 
Fourier space the local extraction involves a parity switch. However, this feature can be trivially removed by defining 
V a = curl curl X a instead of Eq. (26). 

B. rank-2 tensors 

For rank-1 tensors the above considerations are trivial to find and well known. For rank-2 tensors, on the other hand, 
the local decomposition into scalar, vector and tensor modes is not quite so easy. A general projected, symmetric, 
and trace-free rank-2 tensor has the decomposition 

X a b = Sab + V a b + Tab 

= V { aVb)S + V {a Vb } +Tab, (27) 

where the non-local scalar part S a b is curl-free, the vector part V a b is solenoidal, div V = => div div V = (and 
divV a ^ - note the notation: V a V a = => V a V b V ab = V a V b V <Q Vb) = 0), while the tensor part is transverse, 
divT Q = 0. The question is, how do we form local scalar, vector and tensor quantities from X a b and relate them to 
the non-local split given above? That is, what differential operations do we need to do X a b to leave only the scalar, 
vector or tensor part? 

Furthermore, let us assume that X a b obeys a wave equation with source of the form 

C[X ab ] = S ab (28) 

where C contains time derivatives and Laplacians, and any derivative operations which preserve the rank of X ab — 
i.e., curl , dis div or div dis, or combinations thereof. Ideally we would like the local extractions to be differential 
operators which commute with C. To show this we need only show that the Laplacian and curl commute with our 
extractions below: time derivatives commute trivially; dis div commutes if curl and V 2 does by Eq. (17); and div dis 
therefore will by Eq. (18). 

C. scalars 

Clearly, S = div div X is a (covariantly defined) scalar and can only depend on S, so let us define 

Sab = y[X ab ] = V (a V 6) div divX = V (a v fc) V c V d X cd 
= V( a V(,) div div S 

= |V (a V 6> (v 2 +3/v) V 2 S. (29) 
The non-local scalar S is given formally from X a b by 

S = ^(V 2 + 3A") 1 V" 2 divdivX 

= | (V 2 + 3A") 1 V- 2 S. (30) 
This latter relation trivially gives the relation in Fourier space by replacing V 2 M> — k 2 . 
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Defining S ab to preserve the rank of the original tensor allows us to find the wave equation it satisfies easily. First, 
note that any curls in C commute with 5? trivially (producing zero), and note that (non-trivially) 



V 2 (V( a Vf,) div div j = (V< a V b ) div div j V 2 . (31) 

We then see that the locally defined scalar quantity S ab obeys 

C[S ab ] = ,Y[S ab ] = V <a V fc) div div S. (32) 



D. vectors 



We begin by noting that the operator curl div = 2div curl knocks out the scalar and tensor part of X a b, and is 
solenoidal by Eq. (20) (note that it forms a rank-1 tensor from X ab ); thus, V a = curl div X a is a locally defined vector. 
The rank-2 extraction of vector modes may be given by taking the distortion of this operator, 

V ab = f[X ab ] = dis curl d\vX ab 

= e cd(a V b) V c V e X e d 

= dis curl div dis Vab (33) 

= | dis curl (v 2 + 2K)v ab . (34) 

The non-local vector V a may be given formally from X ab by 

V a = 2(-V 4 + 4A 2 ) curl 2 divX a 

= 2 (-V 4 + 4A 2 ) 1 curl V a . (35) 

This last relation enables us to relate the local and non-local vectors in Fourier space as = ^{k 2 + 2K) 1 / 2 {— k 2 + 
2K)V l - k \ / with the same relation for the opposing parity. 

Again, the otherwise superfluous preservation of the rank of X ab in defining the operator y - the extra dis term - 
allows V and C to commute [to prove this for the Laplacian, we use Eqs. (15), (14) and (22); for curl, we use Eqs. (16) 
followed by (10)], giving our wave equation for V ab as 

jO[V ab } = y[S ab ] = dis curl div S ab . (36) 



E. tensors 



The key difficulty here lies in finding the differential operator which when acting on dis V ab = ^ (aYb) produces zero 
by virtue of V a being solenoidal (while obviously leaving the transverse part of X ab ). It is straightforward to verify 
that [curl 2 + ^ dis div — K]A\sV ab = for div V = 0, which follows from Eqs. (10), (12) and (11). We must also take 
an extra curl to remove the scalar part. Therefore, our local tensor extraction may be defined as 



Tab = ^[X ab ] — 



-V 2 + 2A' + 2disdiv 



curl X ab 



+e cde V {a V e V c X b f 



V 2 + 2K ] S7 c X b f + V b) V e V c X 



= -V 2 



2K curlT Qb 



(37) 



Note that curl commutes with the operator in square brackets. It is relatively straightforward to verify that T ab is 
transverse, divT^ = 0, showing that this represents a tensor mode as expected: First note that we may write 



sr 



4curl 



6K 



curl 



(38) 
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which follows from Eq. (17), so that 



V 2 + curl 2 - 2K 



div curl 



i\7div div curl = 0, 



(39) 



which uses Eqs. (16) and (22) on the first line, then Eq. (11) to get to the second, and finally Eq. (9) to show the last 
expression is zero. 

The formal, non-local, TT tensor T ao is given in terms of the original tensor X ao by taking a further curl to obtain 



T ab = (-V 2 + 3Kj 1 (-V 2 + 2Kj 1 #[curl X ab }. 



(40) 



To convert our extraction into Fourier space we define tensor harmonics as V 2 Q^P — — k 2 Q!,T , where we have two 

parities of orthogonal harmonics, (k 2 + 3A') 1/2 QiP = curl Q$ <^ (fc 2 + SK^Q^ = curl Q$ [6]. We therefore 
find 



f« = (k 2 + 3K) 1/2 (k 2 +2K)?( k \ 



(41) 



with the same relation for the opposite parity. 

In order for & to be useful it will have to operate on a wave equation, Eq. (28). It is clear from the form of 2? 
given in Eq. (38) that £ and 2T commute, so that T ao obeys 



C[f ah ] = &[S ab ] = 



-V 2 + 2A' + 2disdiv 



curl5 Qf) 



(42) 



III. VECTOR AND TENSOR MODES INDUCED BY SCALARS 



A useful application of our extraction operators may be illustrated by considering second-order modes in cosmology 
which arc induced by first-order scalar modes. We shall investigate this both in the covariant approach, and briefly in 
the metric based approach. In the former we shall see how the local extraction operators can help in the construction 
of gauge-invariant objects at second-order. 



A. The covariant approach 



Wc shall now consider how to use the relations we have derived in practice, using the 1+3 covariant formulation of 
perturbation theory for illustration [5, 6]. The covariant approach is adapted for perturbation theory because it deals 
with physically defined tensorial objects. Because of this, at a given perturbative order, any tensor which vanishes 
at lower order is automatically gauge- invariant by an extension of the Stewart- Walker Lemma [1, 7] (though they are 
not are not frame-invariant as an observer must be specified). This helps us to use the local SVT operators to form 
gauge-invariant quantities at second-order. 

We introduce an index-free formulation of the covariant approach, as described in the Appendix. In this section 
the spatial derivative operators are the normal ones used in the covariant approach, and do not commute with the 
time derivative. 

Consider a curved FLRW background with some fluid, with an equation of state p = wp. Now perturb this at 
first-order exciting scalar modes only. This can be described by the shear and electric Weyl curvature only: 

er — dis A = -2Ha - E , (43) 
E = -\{l + w)pa - SHE, (44) 



where the acceleration is related to E by 



A = --—^divE. (45) 
(1 + wjp 



The solution to this system determines all other gauge-invariant quantities, such as DiJ = | diver, Dp = 3div_E, 
etc. All these quantities are gauge-invariant by the Stewart- Walker Lemma. The fact that only scalars are allowed 
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is actually signalled by the condition that curler = = curlE, and so div curler = = div curl £7, which all follow 
when the vorticity is zero (at first-order). 4 

If we now allow for second-order perturbations, then er and E are no longer gauge-invariant, and contain a mixture of 
scalars, vectors and tensors, all sourced by the scalars at first-order (that is, the solutions to the first-order equations 
above). However, given our method of a local SVT extraction, we can form objects which are gauge-invariant at 
second-order: the rank-2 tensors ^(er), Y(E), &{cr) and &{E) are all gauge-invariant at second-order by the 
Stewart- Walker Lemma, since they vanish at first-order. It is important to note, however, that these objects may not 
necessarily describe vectors and tensors because er and E do not vanish at first-order, a key requirement in our proofs 
(e.g., 3T{E) may not be divergence-free). 5 A corollary to this is that because "V(ct\ Y(E), 2T{ct) and 3^(E) form a 
set of second-order gauge invariant rank-2 tensors, we can construct pure scalar, vector or tensor modes from them 
by operating further with 5? , Y or & respectively. For example, ^[^(cr)] is a second-order gauge-invariant scalar 
degree of freedom, and f7\y{E)\ is a second-order gauge-invariant tensor mode. These are quite cumbersome to use, 
however, because the objects they are formed from - er and E - are not gauge-invariant at second-order. 

The magnetic Weyl curvature, on the other hand, is a fully gauge-invariant variable at second-order because it 
vanishes at first-order, and so is ideal for local SVT decomposition at second-order. We have 

H = - dis w + curl er (46) 

which tells us that this contains both vectors (signalled by the vorticity) and tensor modes ( curl er is zero for scalars at 
first-order, so gauge-invariantly signals tensor modes at second [8]). Does H contain scalars? No. Consider Eq. (A15), 

div if = - (1 + w) pu) - er x E ; (47) 

when we take the divergence of this, we find, using Eq. (A32) 

div div H = — (1 + w) p div u> + cr • curl E — E ■ curl er (48) 

and this is identically zero at second-order. If we are to consider vector modes, then there are two sources: one is 
from the vorticity, but this isn't a generated vector mode at second-order, since it is straightforward to show that 
d> = (3w ~2)u> [13], on using curl A = 6wHu>. Let us set this to zero. The other vector degree of freedom may be 
found by locally extracting it from if; using the results of Sec. II D, and Eq. (A33) we have: 

curl div H = ^(er • div E - E ■ div er) + er • dis E - E ■ dis er . (49) 

Note that curl div H is a rank-1 tensor describing a pure vector mode which is induced by first-order scalars, as we 
would hope. 

Now let us consider tensors, the tricky part. It is clear that H can represent pure tensor modes only when 
er x E = = uj. More generally, the tensor degrees of freedom are given by a wave equation sourced by scalars. Using 
the identity 

(curlX) = curlX - HcurlX - ^er x div X + er o curl X - cr X dis X + A X X + H A X X (50) 

for E, we find that H satisfies the general wave equation at second-order (i.e., neglecting terms of order H 2 ,crH, 
etc. and keeping uj = 0): 

H-V 2 H + 7HH+ (^6H 2 -2wp+ ( ^pJ H = S, (51) 

where the source <S is given by 
S = - (| + %w) er x dwE - 2 er x dis E + %wer x 6\v E + 3 (-$ + w) E x div er + , - W ' — E x div £(52) 

\b 2 I 2 V 5 J (1+W) p 



4 When vectors and tensors are not set to zero we can use "V and & acting on it or £ to invariantly extract the SVT parts. 

5 This is analogous to the situation at first-order whereby DHor Dp can contain vector degrees of freedom when the vorticity is non-zero, 
despite appearing to be gradients of scalars. Since H and p are non-zero in the background, D is not a pure gradient operator at 
first-order, and so curl Dp <x oj ^ 0. 
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where the div H constraint, Eq. (A15) (which allows a dis div H term to be expressed in terms of first order product) 
is used simplify the first term in the source, which is subsequently expanded using identity Eq. (A32). Now, to get 
the gauge-invariant equation governing the tensor part of this equation it is a straightforward matter of applying the 
tensor extraction operator to this equation. Defining the tensor 



H = a~ 3 ^(H) = [-D 2 + 2a- 2 K + 2dis div] curlH 



(53) 



we have, 



n 



D 2 U 



13 HH 



2 a 



n 



-D' 



2K 
a 



2 dis div 



curl <S, 



(54) 



which locally represents the gravitational wave part of the perturbations. The source terms can be further expanded 
into their invariant symmetric trace-free parts using the identities in the appendix. 



B. The metric approach 



To finish off, let us briefly consider how our approach can be applied in the metric approach to perturbations, in the 
same physical situation considered above where vectors and gravitational waves are sourced by scalar modes. This 
was considered in [9-14]. In the Poisson gauge, using conformal time 77, 



ds 2 = -a 2 [1 + 2$] drj 2 - a 2 Vidx l drj + a 2 [(1 - 2$) 7ij + h l3 ] dx l dx j 



(55) 



where it is assumed that $ is first-order (and there is no anisotropic stress) and Vi and hij are the second-order vector 
and tensor degrees of freedom in the metric. From the ij part of the field equations we have 

2d {l V^ + md {l V 3) + h'l 3 + 2-Hh'ij - V 2 hij = Stj (56) 

where TL = a' /a, a prime denotes d„, and the source is given, for an arbitrary constant equation of state w, by 

16 



-16$%^-)$ -8d (i $3j-)$- 



3(i + w )n : 



-d {i (*' + H®) <%> ($' + W*) 



(57) 



The vector and tensor parts typically have to be separated by invoking an extraction operation which involves a 
non-local projection in Fourier space. For example, to extract the divergence free vector Vi from dnVj}, we can use 
the operator [12], 



2i 



(2tt) 3 



d 3 fc' k'~ 2 / d 3 a;' k H [e m {k')e j {k') + e m {k')e j {k')] e lk '< a 



(58) 



where e and e are vectors orthogonal to fe in Fourier space. Then V m (x) = V%(x,x')duVj)(x'). A similar operation 
isolates the tensor parts in a similar fashion: decompose the object in Fourier space, project out the tensor degrees 
of freedom, and then reconstruct a tensor object in real space. With such operators Eq. (56) can be decomposed into 
two equations for each mode. 

To perform the same operation locally is impossible. However, to use the formalism presented here is straightforward, 
and gives an alternative way to separate Eq. (56) into a vector equation and a tensor one. Let us consider how to 
form the vector equation. Let us first convert Eq. (56) into index-free notation: 



2d\sV' + 4Hd\sV + h" + 2Hti - V 2 h = -16$dis V$ - 8V$ o V$ 



16 



3(l + w)H 2 

Now, if we operate with y on both sides we have an evolution equation for ^(dis V): 



V($' + H$)oV($' + W) 



a 2 r(dis V) 



a 2 y 



-8$dis V$ - 4V$ o V$ 



#' + M)oV($' + H$) 



(59) 



(60) 
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The right hand side can be converted into irreducible form as follows. First note that ^($dis V$) = — Y(V^ o V'P), 
which implies that each term on the rhs is of the form "^(V$ o V\&). Using the identities in the appendix, together 
with the fact that curl V$ = 0. we hnd that: 



dis curl div (V$ o V*) 



5 - 

9 
1 

To 
i 

3 



VV 2 $ x dis W + VV 2 * x dis V$ - V$ x dis VV 2 * - W x dis VV 2 $ 

dis V<1> x div dis V* - dis W x div dis V$ 
dis V<3> x dis dis W - dis V* x dis dis V$ 



(61) 



Note that the last two lines vanish in the case when "J = $. Hence, in the case of dust, or purely growing mode initial 
conditions where oc <E>, these do not contribute to the induced vector mode. In these cases we have, 



-T(disV) 



40 



27(1 + w)a 2 



VV 2 $ x dis V$ - V$ x dis VV 2 $ 



+2VV 2 ($' + H<S>) x dis V($' + "H$) - 2V($' + %<$>) x dis VV 2 ($' + %<$>) 
A similar equation for the tensor part may be found in an analogous fashion. 



(62) 



IV. DISCUSSION 



We have presented a novel way to invoke the scalar, vector and tensor split at the heart of cosmological perturbation 
theory in a local fashion without relying on boundary conditions. Although it appears calculationally cumbersome, 
this is conceptually important in cosmology because the boundary conditions are fundamentally unknowable. We 
have also developed a fully index-free 'vector-calculus'-likc approach for the 1+3 covariant approach which makes 
complicated tensor equations simpler to look at and much easier to manipulate. 

As an application we have examined the case of vectors and tensors which are induced at second-order from first- 
order density perturbations, and have seen that 'f'(H) and S?{H) arc gauge-invariant variables which are a pure 
vector and tensor mode respectively. The covariant approach has a significant advantage over the metric approach 
in this regard: it is much harder to isolate gauge-invariant quantities at second-order starting from a non-covariant 
split of the metric, although we have discussed how to extract local SVT quantities using the same operations. While 
it is mathematically clear what is going on in this example, the physical interpretation of i / (H) and 3?(H) is not 
obvious and deserves further investigation. It is worth noting that in Fourier space curl div H and 3f(H) are related 
to the usual non-local vector and tensor parts of H by k 3 factors only, so the physical interpretation of the local SVT 
variables may not be much different from the usual ones. 

In this example, we have not considered the more involved problem of second-order scalars induced by first-order 
scalars. H contains no scalar modes in this scenario. However, we have mentioned that we can define second-order 
gauge-invariant scalar modes by noting that ^(er), Y(E), ^(cr) and &(E) are rank-2 tensors at second-order which 
are gauge- invariant (since they all vanish at first-order), so these can be used as the basis to form scalar modes at 
second-order (i.e., y[y{cr)], S fi [y(E)], S fi [^{ar)\ and y{3?(E)]), which necessarily vanish at first-order, and are thus 
gauge-invariant . 

It is a more complicated issue to isolate pure second-order SVT modes locally when scalars, vectors and tensors 
are all present at first-order. At a given perturbative order, the local extraction operations only have the SVT 
interpretation when acting on a quantity that vanishes at lower orders (and is thus gauge-invariant at the relevant 
order if it is a covariantly defined object). One way to approach this is to consider something like dis a;: this cannot 
contain tensors at first-order since it is formed from a rank-1 tensor. Therefore, ^(disu;) must be a gauge-invariant 
object at second-order (provided it is non-zero). This is analogous to having vector modes at first-order contained in 
quantities such as Dp which look like scalars but are not in general. Similarly, quantities such as dis Ddiv div er must 
be rank-2 tensors consisting of pure scalars at first-order, and so ^(dis Ddiv diver) and ^(dis Ddiv diver) must 
be gauge-invariant objects at second. In a similar fashion, one can form div d\vY{cr) or div div ^(er), which must 
vanish at first-order and so be gauge-invariant scalars at second-order. In an accompanying paper this is extended 
into a general technique for defining gauge-invariant SVT objects at any order [16]. 
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Appendix A: The 1+3 covariant approach in index-free notation 



Index free notation may be used for any equations which are irreducibly split, and all objects appearing are similarly 
split [15]. In this case we denote a 3-vector V a by V and, more generally, a PSTF tensor X a ...h by X. This notation 
is unambiguous provided the valance of an equation or variable is known. 

We define three products between vectors and PSTF tensors. For this, let V, W be rank-1, X,Y be rank-2 and 
<&, \l/ be rank-3 (rank-3 objects commonly appear as distortions of rank-2 tensors): 



Type 


Rank 


Rank 1 


Rank 2 


Rank 3 


Rank 4 


dot 

product 


V a W a = VW 
X ab Y ab = X Y 
^ abc <f abc = * * 


VX ab = V ■ X 
$ abc X bc = * X 


$ abc V c = * V 






cross 
product 




£ab<y b w c = V X W 

e abc X b d Y cd = XxY 

eabc^ef^ = * X * 


e cd{a V c X d b) = VxX 
e od{a X ce $ d b)e = X x * 


e de(a V d ^ bc) = V x * 
€ de{a X d b Y e c) =XxY 


eef(aX e b <f> f cd) = X x * 


circle 
product 






V (a W b) =VoW 
X (a c Y b)c =XoY 

<f>cd(a i & Cd b) = 


X (ab V c) =XoV 

$(ab d X c)d = * O X 


X (ab Y cd) =XoY 

<$>(ab C V d ) = * ° V 



Note that the notation is unambiguous once we know the rank of the object. For example, $ o \J> can mean a rank 
2, 4 or 6 tensor. However, since this notation is only used on PSTF tensors, once the rank is specified, the product is 
unique. (Indeed, we don't actually need the new symbol o as ■ would do just as well - we have introduced it for extra 
clarity.) Note also that the cross product is anti-symmetric on its arguments. 

In the fully non- linear equations of the covariant formalism we normally use D a rather than V Q , which is only 
really useful on an FLRW background. Let us also define: 

V = aD, (Al) 

curl = a curl, (A2) 

div = adiv , (A3) 

dis = adis . (A4) 

With this notation, all the identities above (and below) which depend on spatial commutation relations, are the 
same, if we replace K h-> — i 3 i? = K/a 2 on an FLRW background. The only difference comes when commuting time 
derivatives. 

The 1+3 evolution and constraint equations are in this notation: 
Evolution equations: 
Rank 0: 

3H - div A = -3H 2 + A ■ A - a ■ a + 2uj ■ u> - \ (p + 3p) + A (A5) 



p+divq = -3H(p + p) - 2 A ■ q - a ■ tv (A6) 

Rank 1: 

w-icurlA = -2Hu> + w • a (A7) 



q+Dp+div7r = -AHq - q ■ a - (p + p) A - A ■ tt - u x q (A8) 

Rank 2: 

<r-disA = -2H<j + A o A - a o er - uj o u - E + ±7r (A9) 



2 



TV - curl H + \ dis q = -\ (p + p) er - 3H [E + ±7r) 

+ 3croE-^croTT-Aoq + 2AxH + ujxE + ^LOX7v (A10) 



H+ curlS - | curl tt = -3HH + 3 cr o H + |w o q - 2 A x E - \q x <r + w x H (All) 
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Constraint equations: 
Rank 0: 



Rank 1: 



Rank 2: 



= cliv u) — A - u) (A12) 

= diver - 2DH+ curlw + 2 A x u + q. (A13) 

= div£;+ idivTr - ^Dp + Hq - ±<x • q - 3u; • H - cr x H + §u; x q (A14) 

= div if + ± curie/ + (p + p)u + 3w • E - |w • tt + cr x £ + ±er x tt (A15) 

= ff + disu; - curler + 2 Aou; (A16) 



In this notation it is actually much easier to manipulate the equations, and write down the various commutation 
relations. For example, 

(div V) = div V + A ■ V - Hdiv V - cr ■ dis V - w • curl V + 2HA ■ V 

A ■ cr -V — A ■ (oj x V) — V ■ div cr + 2V ■ DH - V ■ curl w , (A17) 

is a fully non-linear commutation relation. 

1. Identities 

The following identities are useful for manipulating and expanding derivatives of products. They extend the usual 
rules of vector calculus to PSTF tensors. For example, the first three are normally given as vector calculus identities; 
here we have expanded them further using dis. We use V, W for rank-1 tensors and X,Y for rank-2. We didn't 
quite have the energy to go further but products of rank-3 are actually required for full decomposition of source terms 
in rank-2 wave equations. Note that these are fully non-linear identities as derivatives are never commuted. This 
means that they apply equally well for D as V. 

V(V-W) = v a (v b w b ) 

= i (V div W + Wdw V) - - (V X curl W + W x curl V) + V ■ dis W + W ■ dis V, (A18) 

div (VxW) = V a (e abc V b W c ) 

= W cur\V -V ■ curl W, (A19) 
curl (VxW) = e abc V b (e cde V d W e ) 

= - (Vdiv W - Wdw V) + -(V x curl W - W x curl V) - (V ■ dis W - W ■ dis V) (A20) 
3 2 



dis (VxW) = V (a [e b)cd V c W a 



- (V o curl W-Wo curl V) + V x dis W - W x dis V (A21) 
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div (VoW) = V a (V {a W b) ) 



^ (Fdiv W + Wdw V) - -L (V x curl W + W x curl V) 



-(V- d\sW + W- d\sV) (A22) 



1 



curl (VoW) = ^e cd(a V c (V b) W d + V d W b) ) 

= |(Fo curl W + Wo curl V) - \ (V x dis W + W x dis V) (A23) 

dis (VoW) = V (a (V {b W c)) ) 

= V o d\sW + Wo d\sV (A24) 

div (V ■ X) = V a (V b X ab ) 

= X ■ d\sV + V divX. (A25) 
curl (V ■ X) = e abc V b (V d X c d ) 

= —X curl V + V ■ curl X - X x dis V (A26) 

dis (V ■ X) = V (a V c X b)c 

= -^Xdiv V - \x x curl V + X ■ dis V + ^-V o div X - \v x curl X + V ■ dis X, (A27) 



div (V x X) = V a (e cd<a ^ c X rf fc) ) 



-X • curl V - iv ■ curl X - ^X x dis V + - V x div X, (A28) 
4 2 2 2 

2 fc cd(aV \^| e/ | K ^ 6) 

ixdivF + ixcurl V - ^Fcur! X + ^X o disF - ]-V ■ disX, (A29) 



curl (VxX) = -e cd{a V c (e d V e X f b) +e lem V e X d f 



dis (V xX) = V/ c (e ldel{a V d X e 



2 4 3 2 2 

<c (t\de\(ay d X e b - } ^ 

\v o curl X + X o curl V + V x dis X - X x dis V, (A30) 



V(X-F) = V c (X ab Y ab ) 

= - (div X • F + div Y ■ X) + - (X x curl 1" + F x curl X) + X • dis F + Y ■ dis X, (A31) 
5 3 

div (X x Y) = V Q (e abc X b d Y cd ) 

= X curIF + Y ■ curIX, (A32) 
curl (X x F) = e abc V b (e cde X df Yj) 

= (X • div F - F • div X) + - (X x curl Y - Y x curl X) - X • dis F + F • dis X, (A33) 



dis (X x F) = V (a (e b)cd X^Y d ) 



— (X x div F - F x div X) + - (X o curl F - F o curl X) - X x dis F + F x dis X,(A34) 
10 3 
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div (XoY) = V a (x {a d Y b)d ) 



— (X ■ div Y + Y ■ div X) - — (X x curl Y + Y x curl X) + - (X ■ dis Y + Y ■ dis X) , (A35) 
20 v ' 18 y ' 6 v ' y ' 



curl (XoY] = £cd(a V c Qx de n>e + ^V^e) 



3 11 
-— (X x div Y + Y x div X) + - (X o curl Y + Y o curl X) - - (X x dis Y + Y x dis XA36) 
20 6 2 



dis [XoY) = v (a (^q b y c))d 

= ^- (X o div y + Y o div X) + i (X x curl Y + Y x curl X) X o dis Y + Y o dis X (A37) 

i-U o 
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